Formation of spatial shell structure in the superfluid to Mott insulator transition 
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We report on the direct observation of the transition from a compressible superfluid to an in- 
compressible Mott insulator by recording the in-trap density distribution of a Bosonic quantum gas 
in an optical lattice. Using spatially selective microwave transitions and spin changing collisions, 
we are able to locally modify the spin state of the trapped quantum gas and record the spatial 
distribution of lattice sites with different filling factors. As the system evolves from a superfiuid 
to a Mott insulator, we observe the formation of a distinct shell structure, in good agreement with 
theory. 

PACS numbers: 03.75.Hh 03.75.Lm, 03.75.Mn 
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Ultracold atoms in periodic potentials have proven to 
be highly tunable model systems for investigatingfun- 
damental problems of condensed matter physics 
Among those, the transition from a superfiuid (SF) phase 
to a Mott insulating (MI) phase has attracted much at- 
tention (see P, Q for recent reviews and references) . So 
far, experiments have relied mostly on the time-of-flight 
technique to probe the many-body properties of the gas 
cloud. Unfortunately, one of the most important char- 
acteristics of a MI state, its incompressibility, remains 
hidden in such measurements. In an overall confining 
potential, as is always present for trapped gases in a lat- 
tice potential, incompressibility results in the formation 
of a shell structure. This structure consists of a succes- 
sion of large concentric "Mott plateaus", where the on-site 
occupation is pinned to integer values with suppressed 
fluctuations and abrupt interfaces between two of these 
shells miSHIlllElIIl- In contrast, the conden- 
sate density distribution in the weakly interacting regime 
mirrors the smooth shape of the confining potential 
This dramatic change has so far eluded direct observation 
[i^ . although earlier experiments support the existence 
of a shell structure |Il[llEl- 

Here we report on a direct observation of the in-trap 
density distribution with a high spatial resolution. This 
is achieved by locally modifying the spin state using a 
technique similar to magnetic resonance imaging and to 
the RF addressing technique in j^. Spin changing col- 
lisions ^3 allow us to independently measure the in- 
tegrated density profiles both of the total atom distri- 
bution and of sites with a specific occupation number. 
From such integrated profiles we obtain clear signatures 
of the emergence of Mott plateaus with uniform occupa- 
tion above the MI transition. The measured plateau radii 
agree well with a simple model assuming a fully incom- 
pressible system, zero temperature and zero tunnelling. 

Our experiment can be outlined as follows: After 
preparing the gas in the optical lattice at the desired 
depth, the depth is increased rapidly in order to freeze 
out the spatial distribution in the trap. The atom cloud 
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FIG. 1: Density distribution of a superfluid (a) and a Mott 
insulating state (b) in a lattice with harmonic conflnement. 
Grey solid lines denote the total density profiles, blue (red) 
lines the density profiles from singly (doubly) occupied sites. 
A vertical magnetic field gradient is applied which creates al- 
most horizontal surfaces of equal Zeeman shift over the cloud 
(dashed lines in inset). A slice of atoms can be transferred 
to a different hyperfine state by using microwave radiation 
only resonant on one specific surface (Coloured areas in in- 
sets). Spin changing collisions can then be used to separate 
singly (blue) and doubly occupied sites (red) in that plane 
into different hyperfine states. 



is situated inside a magnetic field with a vertical gra- 
dient and correspondingly almost horizontal surfaces of 
equal field strength as depicted in Fig. ^ A microwave 
pulse on a field-sensitive transition can now selectively 
address atoms from only a sHce of the cloud where the 
transition is shifted into resonance by the field. Atoms 
are transferred to a different hyperfine state within the 
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spectral width of this pulse, corresponding to a region 
with a thickness d of typically a few micrometers. This 
allows the determination of the fraction of atoms in the 
sHce using spin state-selective detection. Repeating the 
experiment for different microwave frequencies yields a 
complete profile of the atom density along the direction 
of the field gradient (Grey curves in Fig. QJ. 

In order to distinguish doubly occupied sites from oth- 
ers, we transfer the atoms of the slice to a hyperfine 
state in which they can undergo spin-changing collisions. 
A controlled spin oscillation then allows atoms on the 
doubly occupied sites to make a transition to another, 
previously unoccupied, spin state (red areas in the inset 
in Fig.^), not affecting sites with occupation numbers 
other than n = 2 (see [18] for details). By subsequently 
measuring the population in each of the states, we obtain 
the spatial distribution of the overall density and of the 
density of doubly-occupied sites independently. 

To discuss the expected shapes and distinct features 
of the resulting integrated density profiles, let us first 
consider the weakly interacting regime of a BEC. In this 
case, the full 3D distribution of the mean atom num- 
ber across the lattice is the Thomas-Fermi distribution 
n-TF — no max(0, 1 — (a;^ + -I- z^)/i?^p), where i?TF is 
the Thomas-Fermi radius and ng the central 3D density. 
The integrated profiles are obtained through integration 
over the x— and y— axes and a sHce of thickness d on the 
2;— axis (see Fig. 1). The integrated distribution of sites 
occupied by a specific number of atoms can be readily 
obtained from the 3D mean atom number distribution, 
by assuming a Poissonian atom number statistics on each 
site before carrying out the integration over the sHce. In 
this limit, it is found that the distributions of both dou- 
bly and singly occupied sites have a smooth shape very 
similar to the profile of the overall density, and with a 
similar width (see Fig.QJi). Doubly occupied sites can in 
fact be found at any place inside the cloud, all the way 
out to the border. 

In contrast, the 3D density distribution of a Mott in- 
sulator state in the low atom number limit for given 
trap parameters is expected to have a spherical shape, 
with constant unity filling up to a radius R. At this 
radius, the density abruptly drops to zero when assum- 
ing negligible tunnehng and zero temperature. The in- 
tegrated density profile of such a sphere with constant 
density is given by an inverted parabola of the form 
j/(z,i?) = SttAj"^ ■(i-max(0, i?^ — z^), where Aiat is the lat- 
tice wavelength. The radius of this parabola R is equiva- 
lent to the radius of the spherical density distribution in 
3D, and the distribution has sharp edges at this distance 
from the center. For higher atom numbers, a two-shell 
system is expected to develop, with an inner spherical 
core of doubly occupied sites and radius R2 and an outer 
shell with occupation number n = 1 of radius Ri. Here, 
the integrated density profile of the n = 2 core is given by 
a parabola h'2{z) = 2i/(z,i?2) (dash dotted curve in Fig. 
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FIG. 2: Integrated in-trap density profiles of the atom cloud 
for different lattice depths and atom numbers: (a) 1.0 x 10^ 
atoms in the superfiuid regime (Vo ~ S-Er), (b) 1.0 x 10^ 
atoms in the Mott regime (Vo = 22E^), (c) 2.0 x 10^ atoms, 
(d) 3.5 X 10"" atoms. The grey data points denote the total 
density distribution and the red points the distribution of 
doubly occupied sites. The blue points show the distribution 
of sites with occupations other than n = 2. The solid lines are 
fits to an integrated Thomas- Fermi distribution in (a), and 
an integrated shell distribution for (b) to (d). The n = 2 
data points are offset vertically for clarity. 



lb). The corresponding n = 1 shell is hollow, and leads 
to an integrated distribution vi^z) = v(z, Ri) — iy{z, R2). 
This is a truncated flat top parabola (dashed curve in 
flgure lb), in stark contrast to what is observed in the 
weakly interacting case. Note that in the MI case also 
the radii of the integrated profile corresponding to singly 
occupied sites and the one corresponding to doubly oc- 
cupied sites are markedly different. 

We prepare our atomic sample by loading a quasipure 
BEC with up to 5 X 10^ ^^Rb atoms prepared in the 
\F = l,mF = —1) state in a QUIC magnetic trap [l9| . 
Combined with a homogeneous offset field, this results in 
an almost isotropic magnetic trap with a frequency ujt ~ 
2tt X 15 Hz for the three axes. Because of gravity, the 
atoms are displaced from the magnetic field minimum, 
leaving them in a gradient field of dB/dz = 3.4G/mm 
which provides the spatial selectivity of the microwave 
pulse. The curvature of the surfaces of equal magnetic 
field due to the lateral confinement of the trap has no 
significant effect on the integrated density profiles within 
our measurement accuracy. The condensate is loaded 
into an optical lattice created by mutually orthogonal 
standing waves with a wavelength of Aiat = 843 nm and 
an average waist of 144 ^m by ramping up the potential 
depth to a final value Vb (typically 22 E^) within 160 ms 
|2(l |. In this configuration the external confinement is 
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provided by a combination of the magnetic trap and the 
Gaussian profile of the lattice beams. For analyzing the 
density distribution, the lattice depth is then rapidly in- 
creased to 40 i?r within 100 ^s in order to freeze out the 
distribution. Here, = is the recoil energy 

of the lattice light, with m the mass of a single atom. 

A microwave 7r-pulse of 150 /is length is then used to 
transfer one slice of atoms from \F — l,mp — —1) to 
|F = 2,TO_F = 0). A second 7r-pulse is applied to bring all 
atoms from \F — 2, mp — 0) to |F = 1, mp = 0) within 
5.6 ^B. The magnetic gradient field is then switched off 
in favor of a homogeneous offset field of 1.2 G. At this 
magnetic field, spin-changing collisions within the F = 1 
manifold are strongly suppressed, even on lattice sites 
which contain more than one atom ■ Using microwave 
dressing 0|, the coUisional spin dynamics can be fully 
controlled and doubly occupied sites selectively trans- 
ferred from the \mF = Q,mp — 0) two-particle state to 
the \mF = —l,mF ~ +1) state. Finally the optical lat- 
tice is switched off and, during a time-of-flight period of 
13 ms, the spin states are separated by a 3.5 ms Stern- 
Gerlach gradient field pulse. The separated atom clouds 
are imaged using standard absorption imaging technique. 

In order to obtain the integrated density profiles from 
the resulting set of images, the atoms corresponding to 
the three Zeeman sublevels are counted by integrating 
over the cloud areas of the absorption images, yielding 
the population numbers Nmp=-i, Nmp=o and Nmp=+i 
for the three m^-states, respectively. Atoms from outside 
the sHce are in the mp = —I Zeeman sublevel, atoms 
from inside the slice which were not in doubly occupied 
sites are in the mp — sublevel. Due to the spectral 
profile of the slicing pulse, the transfer of both atoms of 
a pair has an overall efficiency of a = 0.58 across the 
slice. After the collisional interaction, one atom from 
each transferred pair is in the mp = +1 state and the 
other in the mp = —1 state. The number of atoms among 
those in the slice which originated from doubly occupied 
sites is therefore 1^2 — 2- N„ij,=+i/ a. The total number of 
atoms addressed by the slicing pulse is v^Mce = 
2 • Nmp=+i and iVtotai the overall number of atoms in 
that image. 

Some resulting density profiles can be seen in Figure|2l 
Using the known field gradient and magnetic moment, we 
translate microwave frequency into position. We plot the 
fraction of the atoms in the slice I'siice/A'totai as well as the 
fraction from doubly occupied sites, f2/A'^totai and their 
difference. This difference corresponds to the number 
of atoms from sites with occupation numbers other than 
n = 2 (mostly 1 and 3 for our parameters) . In the weakly 
interacting state (Vq = 3 i^r), the overall profile is fitted 
well by an integrated Thomas-Fermi profile with central 
density and outer radius as fit parameters. The same 
function was fitted to the profile of the doubly occupied 
sites, leading to a very similar radius (see Fig. 

For deeper lattices (Vb = 22E^), the measured inte- 
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FIG. 3: Fitted radii of the Mott shells for different atom 
numbers and a lattice depth of 22 E^. The solid, grey and 
hollow data points denote the radii of the shells with single, 
double and triple occupation, respectively. The solid, dashed 
and dashed dotted lines are calculations of these first three 
Mott shell radii for zero temperature and zero tunnelling. 

grated density distribution is strikingly different from the 
one in the weakly interacting regime. For an overall atom 
number of 2 x 10^ atoms, we observe that the distribution 
of doubly occupied sites is contained in a much smaller 
radius. Moreover, the shape of the number distribution 
is that of an inverted parabola with sharp edges instead 
of the smoother one found for low lattice depths. The 
difference signal between the profiles for the integrated 
overall atom number distribution and the n = 2 distri- 
bution exhibits a truncated parabola profile, as is ex- 
pected for a 3D shell structure with an rt = 1 and n — 2 
shell being present in the trap. In order to quantify the 
observed profiles, we fit parabola shapes according to a 
shell model with two (for small atom numbers) or three 
shells (for the profile corresponding to an atom number 
of 3.5 X 10^ atoms) to our data. The radii of the shells 
as well as the sHce thickness d are fitted independently. 
We observe that with increasing atom number, the size 
of the inner n — 2 Mott shell grows, until for larger atom 
numbers a third shell is formed in the center of the cloud, 
which we observe as an additional parabola on top of the 
truncated parabola of n = 1 sites in Fig. |2li. 

The evolution of the sizes of the Mott shells when in- 
creasing the atom number is shown in Fig Here we 
plot the independently fitted shell radii versus the mean 
atom number of the cloud. The expected radii have been 
calculated for our experimental parameters assuming no 
tunnehng and zero temperature [li|. The distinct growth 
of successive shells with increasing atom number is in 
good agreement with theory, although our experimen- 
tally measured radii seem to be slightly lower than the 
expected values for larger atom numbers. 

Direct observation of the density distribution also 
makes it possible to investigate the evolution of the shell 
structure during the superfiuid to Mott insulator transi- 
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FIG. 4: Evolution of the shell structure through the superfluid 
to Mott-insulator transition for A'^totai = (1-0 ± 0.1) x 10^ 
atoms. Shown is the relative size of the area in which doubly 
occupied lattice sites are found compared to the total cloud 
size. The grey bar marks the lattice depth at which the SF to 
MI transition is expected to occur for a site occupation n = 1. 

tion. For this, we compare the size of the distribution 
of doubly occupied sites to the size of the overall distri- 
bution, and therefore plot the ratio of their respective 
radii Rn=2l RtotaX- Fig. 31 shows the evolution from the 
superfluid state to the Mott insulating state for a flxed 
atom number of A^totai = (1-0 ± 0.1) x 10^. The lattice 
depth was varied between 3 and 40 E^. The corre- 
sponding trap frequencies of the external paraboHc con- 
flnement range from 27rx28 Hz to 27rx80 Hz. To deter- 
mine the radii, we flt a Thomas-Fermi distribution with 
free central density and radius to the respective profiles. 
This distribution is the theoretically predicted one in the 
weakly interacting limit, but was used here as a compro- 
mise for all profiles in order to obtain directly comparable 
values. When crossing the SF-MI transition, we observe 
a drop in the ratio Rn=2/Rtota.\ from « 0.9 below to « 0.6 
above the transition. A theory estimate shows that for 
our trap configuration, between the limits of a weakly in- 
teracting BEC and the limit of a MI, a drop from « 1 to 
~ 0.45 is expected due to the pronounced redistribution 
of the atomic density within the trap. Our data suggests 
that most of the density redistribution happens already 
just before the Mott transition is reached at a lattice 
depth around 13 E^ for our experimental parameters. 

In conclusion, we have observed the density redistribu- 
tion driven by the transition to an incompressible quan- 
tum gas, leading to a shell structure in the in-trap density 
distribution of the Mott insulator state. As temperatures 
on the scale of the on-site interaction energy U strongly 
affect the shape of such a shell structure jlfl l2l] |. this 
method can provide an effective way to measure the tem- 
perature of the system inside the lattice potential. Due 
to the possibility of freezing out the momentary state 
within 100 ^s, the time-evolution of the density distri- 
bution can also be observed. This technique therefore 
could allow for studying the dynamics of excitations in 



the Mott insulator. The local manipulation of the spins 
within the trapped atom cloud with high spatial resolu- 
tion also offers novel possibilities to create and observe 
spin excitations in the lattice system. 
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